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Abstract 

Fusion reactions of neutron-halo nuclei are investigated theoretically with a three-body model. 
The time-dependent wave-packet method is used to solve the three-body Schrodinger equation. 
The halo neutron behaves as a spectator during the Coulomb dissociation process of the projectile. 
The fusion cross sections of 11 Be - 209 Bi and 6 He - 238 TJ are calculated and are compared with 
measurements. Our calculation indicates that the fusion cross section is slightly hindered by the 
presence of weakly bound neutrons. 
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Physics at the dripline is one of current subjects in nuclear physics. Adding neutrons 
to a nucleus as many as possible, we often find that some neutrons are bound very weakly. 
These neutrons form a spatially extended neutron cloud which is called the neutron halo. 
The halo nuclei have been attracting many researchers since its discovery in the secondary 
beam experiment [lj]. The neutron-halo nucleus exhibits many properties different from 
normal nuclei. For instance, it violates the nuclear saturation of density and binding energy. 
Experimentally, these nuclei are mostly studied with reactions using the secondary beam. 

Since the halo nuclei are weakly bound and easily break up, developments of the reaction 
theories capable of describing coupling to the continuum channels are required. For reactions 
in medium and high incident energies, the Glauber and the eikonal theories have been 
successful In these theories, the wave function of the halo neutron in the 

projectile is frozen during the reaction. The success of the eikonal theories rests on the fact 
that the dynamics can be well separated into fast and slow motions. Namely, the relative 
motion of projectile and target is much faster than that of halo neutrons. There is, however, 
no such clear distinction in collisions at low energy. In spite of many experimental and 
theoretical efforts in the last decade, the reaction mechanism of halo nuclei at low incident 
energies is still an open question. 

The fusion of halo nuclei shows an example of the controversial subjects. In early stages, a 
simple and intuitive theoretical argument has been proposed [ill . \v\ : The fusion probability 
was expected to be enhanced because of the spatially extended density of the halo neutron. 
A further enhancement was predicted by the coupling to soft modes inherent to the halo 
structure. Breakup effects were predicted to be small, sustaining a large enhancement of 
fusion cross section especially at sub-barrier energies (l2j, [l3[ . These early studies utilized 
a specific reaction model. We have developed a time-de pen dent wave packet method for 
fully quantum calculations using a three-body model [bl I15L IigL \y\ . Q|. Our results have 
disagreed with the fusion enhancement. In fact, the calculation has even suggested a slight 
suppression of the fusion probability because of the presence of the halo neutron. How- 
ever, recent coupled-channel calculations using a similar three-body model, a g ain , indicate 



a substantial enhancement in the fusion probability at sub-barrier energies 
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Experimental results also show a rather confusing status. In early measurements, an 
enhancement of the fusion cross section at sub-barrier energies was reported for reactions 
using a two-neutron-halo nucleus 6 He . Recently, however, Ref. [2| have reported that 
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the enhancement previously obtained in Ref. |t| is not due to the fusion but to the transfer- 
induced fission. In case of a one-neutron-halo nucleus, 11 Be, the experimental results did not 
present an evident conclusion either jjj. A recent comparison between 10 Be and 11 Be fusion 
cross sections suggests no enhancement for 11 Be [10]. Thus, we may say, at least, that the 
present experimental status indicates no evidence for a fusion enhancement of halo nuclei. 

In this letter, we report calculations of the fusion cross sections in the three-body model. 
We previously reported fusion probabilities of halo nuclei mostl y fo r the case of total angular 
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We now include 



momentum J = which corresponds to the head-on collision 
a full range of impact parameter, for < J < 30, then, for the first time, present results of 
fusion cross sections calculated with the time-dependent wave-packet method. Calculations 
are performed for fusion cross sections of n Be - 209 Bi and 6 He - 238 U. Results are compared 
with experimental data |sl licl]. 

In our three-body model, the projectile is described as a weakly bound two-body system 
of the core and the halo neutron(s). The projectile (P), which is composed of the core 
(C) plus neutron(s) (n), and the target (T) constitute the three-body model. For 11 Be, we 
assume a 10 Be+n structure in which the halo neutron occupies a 2s orbital in the projectile 
ground state. For 6 He, a di-neutron model of a + In is assumed. The time-dependent 
Schrodinger equation for the three-body model is given as 

* s »(H,r,«) = |--v a - 55 v; 

+V„ c (r) + V CT (R CT ) + V nT (r nT )} *(R,r,t), (1) 

where we denote the relative n-C coordinate as r and the relative P-T coordinate as R. The 
reduced masses of n-C and P-T motions are m and //, respectively. 

The Hamiltonian in Eq. (0) includes interaction potentials among constituents. The n-C 
potential, V n c(r), is a real potential. This potential should produce a halo structure of 
the projectile. The core-target potential, Vct(Rct), consists of the nuclear and Coulomb 
potentials. The nuclear potential is complex and its short-ranged imaginary part is inside 
the Coulomb barrier between the core and the target nuclei. The n-T potential, V n T(r nT ), 
is taken to be real. We treat the constituents as spin-less point particles, and ignore the 
non-central forces. The nuclear potentials of V n c, Vct, and V u t are taken to be of a Woods- 
Saxon shape. For n Be - 209 Bi reaction, we take Vo = —53 MeV, r = 1.3 fm, a = 0.7 fm 
for V n c, Vq = —65 MeV, r = 1.16 fm, a = 0.8 fm for the real part of Vct, Vq = —60 
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MeV, tq = 0.6 fm, a = 0.4 fm for the imaginary part of Vct, and Vo = —44.019 MeV, 
r = 1.27 fm, a = 0.67 fm for V nT . For 6 He - 238 U reaction, we adopt V = -69.305 MeV, 
r = 0.667 fm, a = 0.65 fm (V nC ), V = -95 MeV, r = 1.1101 fm, a = 0.5834 fm (real 
Vct), Vo = -10 MeV, r = 0.6 fm, a = 0.7 fm (imaginary V OT ), and V = -93.4 MeV, 
r = 0.986 fm, a = 1.184 fm (V u t)- The V^c provides the 2s bound orbital at —0.6 MeV 
for n Be, and —0.97 MeV for 6 He, which are close to the experimental separation energies of 
one and two neutrons, respectively. The adopted parameter set for V n g is close to that used 
in the studies of the a+ 6 He and p+ 6 He scatterings by Rusek et al. 
gives a large value of B(E2;g.s. — > 2 + ) for 6 He as pointed out in Ref. 

We define the fusion cross section in terms of the loss of flux because of the imaginary 
potential in Vct- The use of the imaginary potential is approximately equivalent to the 
incoming boundary condition inside the barrier. Once the core and the target nuclei are in 
the range of the imaginary Vct inside the Coulomb barrier, the total flux decreases by the 
absorption, no matter whether the neutron is captured by the target or not. Therefore, this 
includes not only complete but also a part of the incomplete fusion in which the charged 
core ( 10 Be for n Be and 4 He for 6 He) and the target fuse while the neutron escapes. 

The wave-packet calculation is performed with the partial wave expansion. Here, we use 
the partial wave expansion in the body-fixed frame 0, 13 Although this is equivalent 
to the one in the space-fixed frame, the body-fixed frame is computationally superior for 
states at finite total angular momentum, J ^ 0. In the body-fixed frame, the channels 
are specified by the projection of J on the body-fixed z-axis, Q, and the magnitude of the 
relative angular momentum conjugate to the angle between r n c and Rpt which we denote 
as I. The coupling between different Q channels is present only for AQ = ±1, given by the 
Coriolis term. This sparse coupling term makes computations much easier in the body-fixed 
channel. Another advantage in the body-fixed frame is a natural truncation scheme set by 
the maximum value of Q. 

We adopt the model space of the radial coordinates, < R < 50 fm and < r < 60 
fm, and discretize them with steps (AR,Ar) = (0.2,1.0) fm (AR,Ar) = (0.3,0.6) fm for 
n Be and and 6 He, respectively. In the case of 11 Be, we also performed a calculation with 
a finer mesh Ar = 0.6 fm and confirmed that the fusion probabilities coincide with the 
present results in a good accuracy. The maximum n-C relative angular momentum is taken 
as Z max = 70 to obtain converged results. If the n-T potential, V n T is switched off, one can 
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FIG. 1: Fusion probability of Be - Bi at J = 0. Calculation with V n T (solid curve) is 
compared with that without V n x potential (dashed), and two-body calculation without a halo 
neutron (dotted). 

use a much smaller cutoff, typically /^ a x = 4. The discrete-variable representation is used 
for the second differential operator |24|. The time evolution of the wave packet is calculated 
employing the fourth order Taylor expansion method. 

The initial wave packet in the incident channel is a spatially localized Gaussian wave 
packet multiplied with an incoming wave for the P-T relative motion while the n-C wave 
function is a bound orbital corresponding to the halo structure. The average P-T distance 
of the wave packet is set as 22 fm and the average incident energy approximately equal to 
the barrier top energy. The time evolution is calculated until the wave packet consists only 
of outgoing waves and the average P-T distance reaches 30 fm. 

A single wave packet solution contains information for a certain range of energy. We 
extract the fusion probability for a fixed incident energy by the energy projection procedure 

0. 

Since the ground state of 11 Be is known to be well described by a single neutron halo 
model, our calculation for n Be is expected to be realistic. On the other hand, a di-neutron 
model for 6 He may be too crude. In order to include two-neutron correlation, we need a 
three-body a + n + n treatment for 6 He. Then, this requires a four-body description for 
the reaction, which is beyond the present computational ability. We consider the present 
di-neutron treatment may provide a qualitative description for 6 He. 

Now, let us show the calculated results. In Figs. [T] and EJ we show calculated fusion 
probabilities at the head-on collision ( J = 0) as a function of the incident energy. In Fig. 
the fusion probability of n Be - 209 Bi (solid curve) is compared with that of 10 Be - 209 Bi 
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FIG. 2: The same as Fig. □ but for 6 He - 238 U. 

(dotted). The latter is a result of a simple two-body calculation with the Vct potential. 
We also show, by a dashed line, a fusion probability of n Be calculated without the neutron- 
target potential, V n T- Figure [T] indicates that the fusion probability is slightly suppressed 
by the presence of the halo neutron. The suppression is observed in both calculations with 
and without V n T- This fact suggests that the suppression originates from the core-target 
Coulomb field which induces the dissociation of the halo neutron. The fusion probability of 
6 He - 238 U system also shows a behavior similar to n Be (Fig.|2j) The suppression of the fusion 
probability is seen to be even stronger than that of n Be. This may be due to a difference 
of the effective charges of the halo neutron: (1/11) (4e) for n Be and (2/6) (2e) for 6 He. The 
Coulomb dissociation effect is more significant for 6 He. There is only a minor difference 
between the solid and dashed curves for both cases, 11 Be and 6 He. We have confirmed that 
the results do not depend on choice of V n T- Therefore, we conclude that reaction mechanism 
associated with V n T does not play a major role in the fusion process of halo nuclei. It may be 
rather surprising that a long-ranged attractive potential of a halo nucleus, which is expected 
in a folding picture of Vct and V n T, does not contribute much to the fusion probability. 
This means that the reaction dynamics are very different from those naturally expected by 
a simple picture. 

In order to understand why, we should look in detail at the time evolution of the wave 

n 

packet. Then, we find the following picture for the dynamics of halo nuclei 16J : When 
the core nucleus is decelerated by the target Coulomb field, the halo neutrons keep their 
incident velocity, leaving the core nucleus. This yields the Coulomb breakup of the projec- 
tile. After the breakup, the fusion takes place between two charged particles, the core and 
the target nuclei. The halo neutrons behave like a spectator. Based on this picture, we 
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FIG. 3: Fusion probability of 11 Be - 209 Bi with different / max . For reference, the fusion probability 
without a halo neutron is plotted by a thin-solid curve. 

propose a mechanism that may explain the reason for the slight suppression of the fusion 
probability. The incident energy of the projectile is initially carried by the core and the 
halo neutron. After the Coulomb breakup, since the dissociated neutrons, as a spectator, 
keep their shared energy, the energy of the core nucleus is smaller than the total incident 
energy of the projectile. This practical decrease in the incident energy makes the fusion 
probability suppressed for a neutron-halo nucleus in comparison with the nucleus without 
halo neutrons. In this spectator picture, one may naively expect that effect of the breakup 
simply leads to an energy shift in curves of the fusion probability. However, our calculation 
indicates that the fusion probability above the barrier is suppressed, while the probability 
below the barrier is almost unchanged. Thus the effect is not that simple and depends on 
the incident e nerg y. 

In Refs. jlfll . I2O] , a strongly enhanced sub-barrier fusion probability was reported for halo 
nuclei by solving a similar three-body problem. This apparently disagrees with our results 
which indicate the suppressed fusion probability. We found that this discrepancy originates 
from a slow convergence of the calculated result with respect to l max [if]]. The treatment 
of the imaginary potential is also different : we include it between the core and the target, 
while it is in the coordinate of the center of mass of the whole projectile relative to the target 
111 Ref. HQ- We examined both treatments of the imaginary potential and found that 
the results are identical. In the coupled-discretized continuum-channels (CDCC) approach 
which is adopted in Refs. 19, 2(j], a small value of / max is applied for the n-C relative angular 
momentum, typically Z max = 4. In Fig. EJ we show the fusion probability for different cutoff 
values of Z max . If we adopt a small / max value, a strong enhancement of the sub-barrier 
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FIG. 4: Fusion cross section of n Be - 209 Bi (solid line) and 10 Be - 209 Bi (dotted line). Measured 
cross sections flfl| are also shown: 11 Be by filled squares and 10 Be by open circles. 

fusion probability is obtained in our three-body calculations. This is consistent to results of 



the CDCC calculation [19|, 20]. However, increasing Z max , the fusion probability gradually 
decreases, and eventually becomes smaller than the values without a halo neutron (dashed 
line). The result of Z max = 30 cannot be distinguished from that of Z max = 50. To get a 
converged result, therefore, one needs to adopt Z max > 30. We confirmed that this slow 
convergence with respect to / max is due to the neutron-target potential, V u t- If we ignore 
V n T and includes only the Coulomb breakup induced by Vct, the convergence is much faster, 
usually / max = 4 is enough. 

Extending the calculations to finite J up to J = 30, we achieve the total fusion cross 
section. In addition to / max , we must also specify a maximum value of Q which we denote 
as fimax- We compare fusion probabilities with f2 max = 0, 5, and 10 for a fixed J, and find 
that the results with f2 max = 5 and 10 are identical. Furthermore, a small difference between 
results of fi max = and fi max = 5 is visible only for energies well above the barrier, E > 42 
MeV. Thus, we here adopt the no-Coriolis approximation (f2 max = 0) to calculate the cross 
sections. 

We show calculated fusion cross sections of 10,11 Be - 209 Bi in Fig. The Vct for 10 Be - 
209 Bi is constructed so that the fusion cross section of this system is reproduced by a simple 
potential picture. Adding a halo neutron, the fusion cross section decreases slightly for an 
entire energy region. This is because the suppression of the fusion probability observed in 
the J = also appears in finite J. 

The measurement indicates that the fusion cross sections of 11 Be do not differ from those 
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FIG. 5: The same as Fig. |I] but for 6 He - 238 U. Experimental data are taken from Ref. a], 
shown by open circles ( 4 He) and by filled circles ( 6 He). Two data points with large error bars at 
Ec.m. < 18 MeV indicate the data for 6 He. 

of 10 Be within the experimental error. At this stage, it is difficult to judge whether or not our 
calculated results of slight suppression in fusion cross section conflict with the measurements. 
In our calculation, a part of the halo neutron proceeds to the forward direction keeping the 
incident velocity after fusion. Therefore, the observation of these forward neutrons in the 
incomplete fusion process may support our picture. In Ref. |H| , the forward neutron emission 
accompanying the fusion is discussed. 

We next consider the fusion cross section of 6 He - 238 U. In Fig. we show the calculated 
fusion cross sections of 6 He and 4 He on 238 U. The experimental data in Ref. [a] indicate, again, 
that there is no difference between 6 He and 4 He fusion cross section within a statistical error. 
Our calculation indicates that the cross section for 6 He is suppressed from the one for 4 He 
at energies above the barrier. The magnitude of suppression in 6 He is larger than that in 
11 Be at high energies. On the other hand, at sub-barrier energies, the calculated fusion 
cross sections of 4 He and 6 He are almost identical to each other. We again need to wait 
for measurements with high statistics to make a definite conclusion whether our calculation 
agrees with the experimental data. 

In summary, we calculated the fusion cross section for nuclei with neutron halo structure. 
The reaction is described as the three-body model, and the three-body Schrodinger equation 
is solved exactly with the time-dependent wave-packet method. We compare calculated 
cross sections with measurements for 11,10 Be - 209 Bi and 6,4 He - 238 U collisions. Recent 
measurements of these systems suggest that there is no evidence for an enhancement in the 
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fusion cross section of halo nuclei. Our calculation indicates that the fusion cross section is 
even suppressed by adding a halo neutron. Further measurements with increased statistics 
are desirable to obtain definite conclusion whether the fusion cross section is increased or 
suppressed by the presence of the halo neutron. 

This work has been supported by the Grant-in- Aid for Scientific Research in Japan (Nos. 
14540369 and 17540231). A part of the numerical calculations have been performed at 
RCNP, Osaka University and at SIPC, University of Tsukuba. 
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